The smallest eigenvalue of Hankel matrices 

Christian Berg* Ryszard Szwarc ^ 
June 24, 2009 

Abstract 

Let TIn = {sn+m),n,m < N denote the Hankel matrix of moments 
of a positive measure with moments of any order. We study the large 
behaviour of the smallest eigenvalue Xn of TIn- It is proved that Aat has 
exponential decay to zero for any measure with compact support. For gen- 
eral determinate moment problems the decay to of Atv can be arbitrarily 
slow or arbitrarily fast. In the indeterminate case, where Atv is known to 
be bounded below, we prove that the limit of the n'th smallest eigenvalue 
of for — > oo tends rapidly to infinity with n. The special case of the 
Stieltjes-Wigert polynomials is discussed. 
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1 Introduction 

Let {sn) be the moment sequence of a positive measure on R with infinite 
support, 

.„^/."*M, „>o. (1) 

By Hamburger's theorem this is equivalent to a real sequence (sn) such that all 
the Hankel matrices 

T^N = {Sn+m)n,m=0^ N = 0,1,... (2) 

*The present work was initiated while the first author was visiting University of Wroclaw 
granted by the HANAP project mentioned under the second author. The first author has been 
supported by grant 272-07-0321 from the Danish Research Council for Nature and Universe. 

^The second author was supported by European Commission Marie Curie Host Fellowship for 
the Transfer of Knowledge "Harmonic Analysis, Nonlinear Analysis and Probability" MTKD- 
CT-2004-013389 and by MNiSW Grant N201 054 32/4285. 



are positive definite. The smallest eigenvalue of ?i at is the positive number 



Xjy = min{{'HNa, a) \ a E C^^^, | |a| | = 1}, (3) 

and clearly Aq > Ai > . . .. The large N behaviour of Aat has been studied in the 
papers [3 El El [IS [221 El] . See also results in [21 [Hj about the behaviour of the 
condition number k{TCn) = ^n/Xn, where A^v denotes the largest eigenvalue of 
Hn- 

Widom and Wilf |22j found the asymptotic behaviour 

A^ ~ AN'/^B'', (4) 

for certain constants A>0,0<-B<lin the case of a measure fi of compact 
support in the Szego class, generalizing results by Szego [H]. In the same paper 
Szego also obtained results about the Hermite and Laguerre case, namely 

A^ ~ AiVl/^5^''^ (5) 

again with certain A, B as above. In all of this paper ~ means that 
o^N I^N ^ 1 as oo. 

Chen and Lawrence [8] found the asymptotic behaviour of Aat for the case 
of /i having the density e"*'' with respect to Lebesgue measure on the interval 
[0,00 [. The result requires (3 > 1/2, and we refer to [8] for the quite involved 
expression. For {3 = \ the asymptotic behaviour is only stated as a conjecture: 

for a certain constant A > 0. 

Chen and Lubinsky [D] found the asymptotic behaviour of Ajv, when /i is a 
generalized (symmetric) exponential weight including e"'^'" with a > 1. 

We recall that the density e"*'' on the half-line is determinate for /? > |, i.e. 
there are no other measures having the moments 

^t-e-''dt = T{^^/f3. (6) 

However, for < /3 < | the density is Stieltjes indeterminate: There are infinitely 
many measures on the half-line with the moments ([6|). The symmetric density 
g-l^l" ig determinate if and only if a > 1. For general information about the 
moment problem see [H [151 US] • 

Berg, Chen and Ismail proved in [3] the general result that the moment se- 
quence dl]) (or the measure /i) is determinate if and only if Aat ^ for ^ oo 
and found the positive lower bound A^v > 1/po in the indeterminate case, where 
Po is given in ([T5l) below. 
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The purpose of the present paper is to prove some general results about the 
behaviour of Aat. 

In section 2 we prove that Xn tends to zero exponentially for any measure n 
of compact support. Theorem 12.31 is a slightly sharpened version, where only the 
boundedness of the coefficients from the three term recurrence relation (171) 
is assumed. We also show that Ajv may tend to zero arbitrarily fast. 

Section 3 is devoted to showing that there exist determinate measures for 
which Atv tends to zero arbitrarily slowly, cf. Theorem 13.61 

In Section 4 we consider the indeterminate case, where Aat is bounded below 
by a positive constant. We prove that the ri'th smallest eigenvalue \N,n of ([2]) 
{n < N) has a lower bound Aoo,n = lim^v^oo Aat „, which tends rapidly to infinity 
with n, cf. Theorem 14.41 To describe our results in detail we need some more 
notation. 

We let (Pn) denote the sequence of orthonormal polynomials with respect to fi, 
uniquely determined by the requirements that P„ is a polynomial of degree n with 
positive leading coefficient and the orthonormality condition / PnPmdfi = 6nm- 

The orthonormal polynomials satisfy the following three-term recurrence re- 
lation 

xPn{x) = hnPn+l{x) + a„P„(x) + 6„_iP„_i(x), (7) 

where &„ > and a„ G M. 

We need the coefficients of the orthonormal polynomials (P„) with respect to 

n 

Pn{x) = Y,KnX\ (8) 
k=0 

and consider the infinite upper triangular matrix 

B = {bk,n), h,n = 0, k> n. (9) 

Let Bjy denote the (A^ + 1) x (A^ + l)-matrix obtained from B by assuming 
k,n < N and let A^^^ = BnB]^. Defining the kernel polynomial 

N N I N 

Kn{z,w) =Y,Pn{z)Pn{w) = \ Yl ^^'^^^'^ | -^'^^ ^^^^ 

we see that A^^^ = {O'fk') is the (A^ + 1) x (A^ + l)-matrix of coefficients to z^w'' 
in Kp^{z,w). The following result going back to A.C. Aitken, cf. Collar [12], has 
been rediscovered several times, see P| IT7]. 

Theorem 1.1. 

A^''^=n],\ 
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For completeness we give the simple proof of Theorem ll.lt 
For < < we have by the reproducing property 



x^Kn{x, y) dfi{x) = y^. 



On the other hand we have 



N N 



j=o e=o 
and therefore 

AT 

e=o 

The following Lemma is also very simple. The identity matrix is denoted 
Lemma 1.2. As infinite matrices we have 

B{B*n) = {B*n)B = B*(nB) = i, 

and B*T-C is an upper triangular matrix. 

Proof. The matrix products B*T-C and HB are well-defined because B is upper 
triangular, and we get 



{B*H)j,k = ^bn,jSn+k = / Pj{x)x^ d^{:, 

n.=0 



which is clearly for j > k, so B*H is also upper triangular. Therefore, B{B*1-L) 
is well-defined and upper triangular. For I < k we finally get 

k j k / k \ 

{B{B*n))l^k = ^ bi,j ^ K,jSn+k = X] I X] bijhnj j Sn+k = 8l,k 
j=0 n=0 n=0 \j=0 J 

by Theorem 11.11 with N = k. 

The relation {B*H)B = B*(T-CB) = / is an easy consequence of the orthogo- 
nality of (P„) with respect to /i. □ 

We also consider the infinite matrix 

/C = (k,,,), = ^ [ P,(e**)Pfc(e-^*) dt. (11) 

Jo 
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It is a classical fact that the indeterminate case occurs if and only if 



J2\Pn{z)\^ <00 (12) 
n=0 

for all z E C It suffices that ( JT2l) holds for just one point £ C \ M, and in this 
case the convergence of f|T2|) is uniform on compact subsets of the complex plane. 

In the indeterminate case we can let ^ oo in ffTOj) leading to the entire 
function of two complex variables 

oo oo 

K{z, w) = J2 Pn{z)Pn{w) = ^ aj^k^'w^ (13) 

ra=0 j,k=0 

and we collect the coefficients of the power series as the symmetric matrix 

A=ia^,k)- (14) 

In Proposition 14.21 we prove that the matrices A, B, K, are of trace class in the 
indeterminate case and 

tr {A) = tr (/C) = po, 

where po is given by 



po = — / K{e'\e-'')dt = — V|P,(e^*)|'dt<oo. (15) 
^vr Jo ^TT Jo 



In the indeterminate case the infinite Hankel matrix Ti = {sn+m) does not cor- 
respond to an operator on i'^ defined on span{5„|n > 0}. In fact, by Carleman's 

theorem we necessarily have Xl^^o ^2^^"^"^ < ^5 hence S2n > 1 for n sufficiently 
large, and therefore 



oo 



Sn^^ = oo for all n. 

m=0 

It is likely that Theorem 11.11 extends to the indeterminate case in the sense 
that AH = TCA = I, where the infinite series ak,isi+j defining AH and HA 
are absolutely convergent. We have not been able to prove this general statement, 
but it holds for the Stieltjes-Wigert case which is treated in Section 5. 

The Stieltjes-Wigert polynomials P„(x;g) are defined in (l50l) . They are or- 
thogonal with respect to a log-normal distribution, known to be indeterminate, 
and the corresponding moment sequence is s„ = Z^. It is known that the 

modified moment sequence (s„) given by s„ = s„ for n > 1 and 



5o = so-(X;Pn(0;g)^)-^ 

n=0 
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is determinate, and the corresponding measure /i is discrete given by 

/i = ^cj^, (16) 

where X is the zero set of the reproducing kernel K{0, z) defined in ( IT3ll and 

c,= (f2Pk{x;q)A , xeX. (17) 
\k=o / 

The Hankel matrices H = (sj+fc) and H = (sj+fc) agree except for the upper 
left corner. In Theorem 15.21 we prove that the smallest eigenvalue Atv of the 
Hankel matrix Ti^ tends to zero exponentially (while A at is bounded below). We 
do it by determining the corresponding orthonormal polynomials Pn{x]q), see 
Theorem 15.31 

2 Fast decay 

We start by proving a lemma which is essentially contained in jH §2]. 
Lemma 2.1. For each G C with \zo\ < 1 we have 

N - -1 



A7V< [{l-\zo\')J2\Pn{Zo) 



n=0 



Proof. For any a G C^+^, a 7^ we have by (j3j) 

{Ti-Nd, a) 



Xn < 



12 



1 1*^! , 

This means that for any non-zero polynomial 

N N 



p{x) = ^ akx'' = ^ CnPn 

we have 



k=0 n=0 



f\p\ dfl 

2^ Jo \p{e''Wdt 



Moreover, by Cauchy's integral formula 

±1 

2tt Jo e** - zq 
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hence by Cauchy-Schwarz's inequality 

|pWI^<^f |p(.«)P*^f (20) 

and the last integral equals (1 — |zoP)~^ by a well-known property of the Poisson 
kernel. Combining (fT9l) and (!20l) for the polynomial 



N 

Pix) = ^PniZo)Pn(X) 
n=0 



leads to 



(i-kor)lm)r \ ^ 



□ 



Remark 2.2. It follows immediately from Lemma [2.11 that if A at > c > for all 

A^, then 



n=0 



for all zq with |2;o| < 1, hence is indeterminate. 

The following theorem proves that Xn tends to zero exponentially in the sense 
that there is an estimate of the form 

\n < AB^, A>0,0 < B <1, (21) 

whenever the measure in ([T]) has compact support. 



Theorem 2.3. Assume that the sequence from (JTj) is bounded with b : = 
limsup6„. Then 

hmsupA^ <Y^^- 

Remark 2.4. Notice that the condition limsup6„ < oo implies that = 
oo, so by Carleman's theorem the moment problem is determinate, cf. [H p. 24]. 
We also recall the fact that fi has compact support if and only if (a„), from 
([7]) are bounded sequences. 

Proof. Taking zq = ai, where < a < 1, we obtain from Lemma [2?T] 

Xn < f (1 - a') l^-MI' j < ((1 - cy'WN-iiat)\' + \P^im)\'])-' . 
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n=0 



Since the distance from the point ai to the support of the orthogonahty measure 
is at least a, we obtain by [201 Remark 2, p. 148] 



hm sup A 



l/N 
N 



< 



26^ 



As a is an arbitrary number less than 1 we get 

252 



lim sup A 



l/N 
N 



< 



1 + 262" 



□ 



Theorem 2.5. For any decreasing sequence (r„) of positive numbers with tq = 1 
and limr„ = 0, there exist determinate probability measures fi for which Xn ^ tn 
for all N. 

Proof. We will construct symmetric probability measures /i with the desired prop- 
erty. Let 

xPn{x) = hnPn+l{x) + 6„_lP„„l(x) (22) 

be the three-term recurrence relation for the orthonormal polynomials associated 
with a symmetric /x. We shall choose 6„ > 0,n > such that Aat < tat for all 

> 0. We always have Aq = tq = 1 because /i is a probability measure. Since 
si = we know that Ai = min(l, S2), S2 = h^, so we can choose < 60 < 1 such 
that Ai = Ti. 

By Lemma 12.11 with = we get 



N 



-1 



A^< \ Y.\Pnm 



vn=0 



\2Af+l 



< A2Ar < 



^2n(0) = (-1)^ 



and in particular 
By ((22]) we have 
and defining 
we get 

A2Af+i < A2Ar < rlr^ 
and we will choose r^, A; > 1, such that 



^2^(0)' 

^6o&2 • • •&2n-2 ^ , 

Olh ■ ■ ■ hn-l 



(23) 



J2k-l 



&2fc-2 ' 
2„2 



k>l 



N>1, 



rlrl...rl <T2N+i, N >l. 
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First choose < ri < and when ri, . . . , tat-i have been chosen, we choose 



< rjv < min 1 



ri . . . tn-i 



It is clear that the sequence (r^) can be chosen such that rfc — > 0. We next define 
hi = Tibo and we finally have an infinity of choices of 62^-1, ^2A;-2 > to satisfy 
rk = hk-i/hk-2-, k>2. 

If (rfe) converges to zero, the decay of A„ is faster than exponential. Clearly 
the corresponding moment problem is determinate since 

\P2nm = {rir2...rn)-' > 1. 

In particular, the unique measure /i solving the moment problem carries no mass 
at 0. □ 

After having chosen the numbers we have several possibilities for selecting 
the coefficients We will discuss three such choices. 

Example 1. For A; > 2 let &2fc-2 = 1 and &2fc-i = and assume that 0. 
Then the corresponding Jacobi matrix J is bounded and it acts on i"^ by 

Let us compute the square of J. We have 

( J^X)„ = bnbn+lXn+2 + + bl)Xn + 6n-2&n-ia^n-2- 

By the choice of we get and + 1. Therefore the operator 

is of the form = I + K, where K is a compact operator. Hence its spectrum 
consists of a sequence of positive numbers converging to 1. Thus the spectrum 
of J is of the form a{J) — {±tn}, where tn is a sequence of positive numbers 
converging to 1, so the measure /i is discrete with bounded support. 

Example 2. Let 62^-2 = rj^^ and &2fc-i = 1 and assume — > 0. Then the 
corresponding Jacobi matrix J is unbounded. By the recurrence relation we have 

x'^P2n{x) = b2„b2„+lP2n+2{^) + i^ln-l + bln)P2n{x) + &2n-2&2n-l-f2n-2(a:;) ■ (24) 

Then Qn{y) = P2n{-\/y) is a polynomial of degree n satisfying 

yQn{y) = r-liQn+i{y) + (i + ^n+i)<5n(y) + r-^Qn-i{y). 

Letting i?„ = and A„ = (1 + t'I^) we get 

d2 2 
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so by Chihara's Theorem (see flUi Th. 8] and [211 Theorem 2.6]) we see that the 
orthogonahty measure u for Qn{y) is discrete. However, v is the image measure 
of the symmetric measure /i under the mapping a; — ^ x^, so also /i is discrete with 
unbounded support. 

Example 3. Let 62^-2 = and 62^-1 = V • With Qn{y) = P2n{y/y) as in 

Example 2 we get from fl2^ 

l/Qn(Z/) = Qn+liy) + anQniy) + <5n-l(Z/) 

where a„ = r„ + l/r„_|_i. If — > we see again that /x is discrete with unbounded 
support. 

3 Slow decay 

The goal of this section is to prove that there exist moment sequences such 
that the corresponding sequence {\n) from ([3]) tends to arbitrarily slowly. This 
is proved in Theorem 13.61 

Consider a symmetric probability measure ^ on the real line with moments 
of any order and infinite support. The corresponding orthonormal polynomials 
{Pn) satisfy a symmetric recurrence relation (122|) . where 6„ > for n > 0. For 
simplicity we assume that the second moment of /x is 1, i.e. S2 = = 1 and 
hence Aq = Ai = 1. This can always be achieved by replacing by dfj,{ax) 

for suitable a > 0. Note that Pq = 1, Pi{x) = x in this case. 

Lemma 3.1. Let (Pn) denote the orthonormal polynomials satisfying ( 1221) with 
bo = 1. The sequence 

Un = \Pnii)\, n>0 (25) 

satisfies 

1 ^n— 1 / \ 

Un+1 = —Un + — Un~l, U > 1, (26) 

On On 

with Uq = Ui = 1. Moreover, for n >0 

\Pn{z)\ < Un, \z\ < 1. 

Proof Let fc„ = bn,n denote the (positive) leading coefficient of P„ and let 
Xi,X2 ■ . . ,Xn denote the positive zeros of P2n- Then 

P2n{x) = k2n{x^ " x\){x^ - xl) . . . {x^ - X^), 

hence 

U2n = (-l)"P2n(0 > 0- (27) 

Similarly, let yi, 2/2 • • • , Vn denote the positive zeros of P2n+i- Then 

P2n+l{x) = k2n+lX{x^ - yl){x^ - yj) ■ ■ ■ (x"^ - Z/^), 



10 



hence 



Combining ([22]), (|27|) and ([28]) gives ([21 
By ([22) we get for \z\ < 1 

1 



W2n+l = (-ir+'^P2n+l(0>0. 



(28) 



\Pn+l{z)\ <-\Pn{z 



bn-1 



Therefore, (l26l) can be used to show by induction that |-Pn(-2)| < Un- 



□ 



Proposition 3.2. Assume that the coefficients (6„) from ( 122]) satisfy bo = 1 and 
bn-i + 1 < for n > 1 and let Un = \Pn{i)\- Then 



max(M2„, U2n+i) < ]^max 
Proof. Since 



k=l 



1 + &2fc-2 1 + &2fc-l 
^2fc-l ' b2k 



n > 0. 



1 + &fc-i 



< 1, 



«fe+i < max(Mfc_i,Mfc), > 1. 
Mfe < max(Mfc_i,Mfc), 
max(Mfc,Mfc+i) < max(Mfc_i,Mfc), k>l. 



we get from 
We have clearly 
thus 

This implies by 

Mn+1 < 7 max(M„_i,u„) < max(u„_2, 

On On 

and replacing n by — 1 in the first inequality 

, 1 + bn-2 I . 

Un < —r max(M„_2,Mn-l)- 

On-1 

Combining the last two inequalities gives 

■l + 6n-2 l+^n-l 



max^Mn, Un+i) < max 



bn-l 



max(M„_2,Mn-i), n>2, 



which implies the conclusion because Uq = Ui = 1. 



□ 
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Lemma 3.3. Let [hn) and (un) be as in Proposition \3.S[ Then the sequence of 
eigenvalues {\n) from ([3]) satisfies 



\k=0 J 



Proof. By [H (1.12)] we have 



2n N ^1 







k=0 

The conclusion follows now by Lemma [3. 11 which shows that |Pfc(e**)| < Uk- □ 



Using the assumption of Proposition 13. 2[ we adopt the notation 

l-r]k= max — , , k > 1. (29) 

V 02^-1 02k J 

Proposition 3.4. Let (6„) and (un) be as in Proposition \3.2[ Then the sequence 
of eigenvalues (Xn) from ([3]) satisfies 



(N k \ "1 

k=l 1=1 J 



Proof. By Lemma 13.31 and the fact that mq = Mi = 1 we have 

(N ^ 
k=l 

Proposition 13.21 states that 

k 

max(M2fc,M2fc+i) < JJ(1 - Vi)- 
1=1 

These two inequalities give the conclusion. □ 
Lemma 3.5. Let (6„) be as in Proposition \3.B and define ^„ by 

bn-i + 1 

7 = 1 - 4n, n>l. 

On 

Then 

n r n— 1 k 

bn = l[{l-Ck)-' 2 + J2ll^^-^i) , n>l. (30) 

k=l L k=l 1=1 
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Proof. We have 

bn = {l- in)-\l + bn-l) = (1 - (l + (1 " Cn-l)"'!! + &n-2)) = ■ ■ ■ , 

and after n steps the formula ends using bo + 1 — 2. □ 

Theorem 3.6. Let (t„) 6e a decreasing sequence of positive numbers satisfying 
Tn — > and To < 1. T/ien ^/lere exists a determinate symmetric probability 
measure /j, onM. for which \n > ^tn for all N. 

In other words, the eigenvalues Xn can decay arbitrarily slowly. 

The proof depends on the following 

Lemma 3.7. Let (e„) be an increasing sequence of positive numbers such that 
Co > 1 and lirne^ = oo. There exists a strictly increasing concave sequence {dn) 
such that do — 1, dn < e„ for all n and limd„ = oo. 

Proof. Define a function f{x) on [0, oo) by /(O) = Cq and f{x) = e^j for n — 1 < 
X < n, for n > 1. This function is left continuous. The discontinuity points 
in ]0, oo[ are denoted by for a strictly increasing subsequence of natural 
numbers. Consider the sequence of points in the plane given by Ao=(0, 1) 
and Ak = {uk, e„^) for > 1. If we connect every two consecutive points Ak and 
Ak+i by the line segment we will obtain a graph of a strictly increasing piecewise 
linear function g{x) such that g{x) < f{x). Moreover g{x) tends to infinity at 
infinity. We are going to construct the graph of a concave function h{x) such that 
h{x) < g{x), h{0) — 1 and h{x) — > oo as x — > oo. Once it is done the sequence 
dn = h(n) satisfies the conclusion of the lemma. We will construct the graph of 
h{x) by tracing the graph T of g{x). The points of T where the slope changes will 
be called nodes. 

We start at the point (0, 1) and draw a graph of the function h{x). We go 
along the first line segment of F until we reach the first node. Then we inspect the 
slope of the next line segment of F. If it is smaller than the slope of the previous 
segment we continue along F until we reach the next node. Otherwise we do 
not change slope and continue drawing the straight line (below F). In this case 
two possibilities may occur. The line does not hit F. Then the graph of h{x) is 
constructed. Otherwise the line hits F. Then two cases are considered. If the line 
hits a node of F, then we follow the procedure described above for the first node. 
If the line hits an interior point of a segment 7 of F, then we continue along the 
segment 7 until we reach the next node, where we follow the procedure described 
for the first node. We point out that the slope of the segment 7 is necessarily 
strictly smaller than the slope of the straight line followed before hitting 7. 

In this way a graph of h{x) with the required properties is constructed. Ob- 
serve that if the graph of h{x) has infinitely many points in common with F, 
then clearly h{x) — > 00 as x — > 00. But if there are only finitely many points 



13 



in common with F, then h{x) is eventually linear with a positive slope, hence 
h{x) — s> oo as X ^ oo. 

□ 

Proof of Theorem \3.6[ Defining e„ = l/r„, there exists by Lemma [321 a concave, 
strictly increasing sequence with do = 1 and limdn = oo and such that dn < 
Cn- Moreover, we may assume that (i„ < n + 1 by replacing dn by min{dn, n + 1). 
In this way we may also assume that d2 < 3. This implies that there exists a 
decreasing sequence of positive numbers Ck,k > 1 such that Ci < 1 and 



2n 



1 + 2 J]cfc. 



k=l 



In fact, we define 



Ci = {d2 - l)/2, C„ = (C?2n - rf2n-2)/2, ^ > 2, 



SO (c„) is decreasing because (i2n is concave. 
Let the sequence 7]^ be defined by 

1 - ^1 = v/c7, l-r]k = 



Ck-l 



k> 2. 



(31) 



n k 



(32) 



Then r]k > and 

d2n = l + 2^J](l-r//)l 

A:=l /=1 

Define the sequence by 

6fc-i = 6fe = ^ > 1- 

Inspired by formula fl5Ul) we finally define a positive sequence (6„) by 6o = 1 and 



k=l 



n— 1 A; 



fc=i 1=1 



n> 1. 



Then we get for n > 1 



J2n 



< 



k=l 

2 + d 
2 — — 

d2n — d 



Vk 



1-2 



2n-2 



< 2 



2n-2 



n— 1 k 

2 + C;2n 
d2n — d2n~2 ' 



^0 
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where we used formula (132|) . This gives 



1 ^ d2n ~~ d2n-2 
2(2 + d2n) 



and since d2n tends to infinity we get 



oo ^ 



CXD. 



In fact, assuming the contrary we get 



oo > 



oo oo 



don — d' 



2n ~ W2n-2 



^ h2n ^ 2(2 + d2r. 
n=l n=l 

SO there exists G N such that for all p G N 



N+p 



don — d' 



2n ~ "2n-2 



2 " 2 + dsn 



> 



E 

n=7V+l 



(ion — d 



2n — "2n-2 _ U'2N+2p " ^2N 



do 



do. 



2 + d; 



2N+2p 



2 + d 



2N+2p 



(33) 



but the right-hand side converges to 1 for p — >■ oo, which is a contradiction. 

The positive sequence (6„) defines a system of orthonormal polynomials via 
fl22|) . The corresponding symmetric probability measure is determinate by Car- 
leman's theorem because of ( 133|) . Moreover, by Proposition 13.41 and formula (l32l) 

we get 

1111 ^ 

— > = — T27V ^ ~T2Ar+l- U 

2d2N 2e2N 2 2 



^2N > ^2N+1 



4 The indeterminate case 

Let (sn) be the moment sequence ([1]). The inequality 



N 



N 



^ s„+ma„am > |afc|^ aGC 



TV+l 



n,m=0 



fe=0 



can be rewritten 



If we write 



N ^ N 

fc=0 fc=0 
AT AT 

k=0 n=0 



(34) 
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and use (j8]), then (IMI) takes the form 



AT N 



N 



n=k 



n=0 k=0 

This immediately gives the following result: 

Lemma 4.1. The eigenvalues Xn are bounded below by a constant c > if and 
only if the upper triangular matrix B = {hk,n) given by ([9]) corresponds to a 
bounded operator on of norm < 1/ \fc. 



Recalling that the indeterminate case was characterized in by Aa? being 
bounded below by a positive constant, we see that the indeterminate case is 
characterized by the boundedness of the operator B. For a characterization of 
the lower boundedness of A^r in a more general setting see [7j. As noticed in [U 
Remark, p. 72], the indeterminacy is also equivalent to the boundedness of the 
matrix /C, cf. ( ITTj) . which is automatically in trace class if it is bounded. 

Concerning the matrices A, /C, B, given by ( fT4l) . (fTTl) .(!9l) respectively, we have: 



Proposition 4.2. Assume that fi is indeterminate. Then the following matrix 
equations hold 

(i) 1C = B*B, 

(ii) A = BB*. 

A, B, K are of trace class and 

tr{A) = ir(/C) =po, 



where po is defined in (fT5l) . 
Furthermore, the sequence 

/ oo \ 1/2 

Cfc = a/o^ = I ^ l&fc.nM , (35) 
\n=k / 

satisfies 

limA;^ = 0, (36) 

fc— ►oo 

and the matrix A = {aj,k) has the following property 

oo 

J2 Wj,k\' < (37) 

j,k=0 

for any e > 0. 
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Proof. From ([8]) we have 

bk,n = ^ 1 1_ Pn{z)z-^'+''> dz = r-'^ j^^ P„(re^*)e-^^* dt. (38) 
Consider r = 1. Then, by Parseval's identity we have 

N n . I.27T N 

EEi^Mp = i/ (39) 



2vr 

n=0 k=0 ^ " n=0 

Therefore, in the indeterminate case the matrix B is Hilbert-Schmidt with Hilbert- 
Schmidt norm p^j"^ , cf. (fTSl) . Hence both B*B and BB* are of trace class with 
trace po- Formula (i) of Proposition 14. 21 is an immediate consequence of Parseval's 
identity. 

We know that Kn{z, w) defined in (ITUl) converges to K{z, w), locally uniformly 
in C^, hence 

N 

4? = E ^j^rih^n aj,k (40) 
?i=max{j,fc) 

for each pair (j, k). The series 

oo oo 
n=max(j,fc) n=0 

is absolutely convergent for each pair (j, k) because B is Hilbert-Schmidt, so (HUll 
implies (ii). 
Defining 

Cfc = ||i3*4||= JZI&mIM , (41) 

\n=k / 

where Sk, k = 0,1, . . . denotes the standard orthonormal basis in i"^, we have the 
following estimate for r > 1 using the Cauchy-Schwarz inequality 



9 oo oo 



'k,n 



^fc=0 / fc=0 fc=0 fc=0 n=k 

However, by (l38l) and by Parseval's identity we have 

oo oo oo n oo _ „27r 

E^''EiM^ = EE^''iM^ = E^/ (42) 

A;=0 n=k n=0 fc=0 n=0 ^ 

Let now 

P{z)={Y^\P^{z)\'\ , zee. (43) 

\n=0 / 
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We finally get 



E-'£:74=T(^ri^(-")l=*)"'<-. 



k=0 

but since 



{m5k,5k)<\\\B*\Sk\\ = \\B*5k\\ 

this shows that \B*\ and hence B is of trace class. 

For a given e > we have P{z) < Cge"^'^' by a theorem of M. Riesz, cf. [H 
Th. 2.4.3], hence by (ST]) and (|42|) 



k=0 



For r = k/e we get in particular 



k 



hence 



4 < cy', 



lim sup A; ^/c^ < ee, 

fe— >oo 

which shows fl36|) . 

Using |aj,fc| < c^Cfc, it is enough to prove that J^'kLo^l < cxo for < e < 1, 
which is weaker than fl36l). □ 



For a sequence a = G £^ we consider the function 

oo oo 
F^{Z) =Y,OlnPn{z) = Y.PnZ'', (44) 
n=0 n=0 

which is an entire function of minimal exponential type because 

\Fa{z)\ < \\a\\P{z), 

where P{z) is given by (|43|) . The following result is a straightforward consequence 

of dBD. 

Proposition 4.3. The sequence of coefficients (3 = of the power series of 
Fa belongs to i'^ and is given hy (3 = Ba. The operator B : i'^ i'^ is one-to-one 
with dense range 



For a compact operator T on 1? we denote by an{T),n = 0, 1, . . . the singular 
values of T in decreasing order, i.e. 

(Jn{T) = min max (45) 

yc£2,dim y=n \\v\\=l, v±V 
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Theorem 4.4. Assume that fi is indeterminate. Let 

Xn = Atv.o < Xn^i < . . . < Xn,n 
denote the + 1 eigenvalues ofTi^ and let 

Aoo,n = lim AAr„. 

N^oo 

For < n < N we have 

1 1 



an{A) = an{B*f > > (46) 



and 



2 



lim yajA) = 0, lim ^ = oo. (47) 

n— >oo n^oo 77,^ 



Proof. By f H5|) we get 



o'n('B*) < max ll-B*?^!!. 

||d||=1, v±So,...,5„-i 



Let n„ denote the projection onto {60, ... , 6n-i}^. Thus by (jH 

a„(i3*) < ||i3*n„|| < ($^c 

\k=n 

On the other hand, for r > 1 we have 

^ ^ (M)V^ g(r) 

k=n k=n 

where 

00 

5(r) := ^{k\ckfr^ < 00 

fc=0 

because of (l36ll and ~ fc/e, which holds by Stirling's formula. Therefore 



.*^2. S{r] 



(n!)2r"' 

and since (JniB*) = we get 



a„(^)=ajSS*)<-^, r>l, (48) 



which proves the first assertion of WJ\ . 
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Let Pr^v denote the projection in onto span{(5o, . . . , 5^}- We then have 



{BPrN){BPrNY = BPr^B* < BB*, 

and therefore for n < N 

an{BB*) > an{{BPrr,){BPrr,y) = cXniB^B^) = an{n^'), 

where the last equahty follows by Theorem ll.li The matrix 7i^^ is positive 
definite, so its singular values are the eigenvalues which are the reciprocals of 
the eigenvalues of TCn, i-e. cr„(H^^) = l/X^^n- This gives (H6l) and the second 
assertion in ( H71) follows. □ 

Theorem 4.5. The trace class operator A : i"^ i"^ is positive with spectrum 

a{A) = {0}U{X^\\n = 0,l,...}. 

Proof. We will consider A^^'^ = i^'fl^) and i^Ar as finite rank operators on i'^ by 
adding zero rows and columns. Clearly, Bn tends to B in the Hilbert-Schmidt 
norm, and therefore A^^^ = B^B]^ tends to .A = BB* in the trace norm. 

The result now follows since the spectrum of A^^'^ consists of the numbers 
A-|^, n = 0, 1, . . . , iV, by Theorem O □ 



5 The Stieltjes-Wigert polynomials 

For < g < 1 we consider the moment sequence s„ = given by 



/27rlog(l/g) Jo ''''^ V 21og(l/g)7 



4^^^ dx. 



(49) 



We call it the Stieltjes-Wigert moment sequence because Stieltjes proved that it 
is indeterminate (he considered the special value q = ^) and Wigert [23j found 
the corresponding orthonormal polynomials 



P„(x;g) = (-ir 



nil n 
q2^4 



\/{<l\<l)n 



fc=0 



(50) 



Here we have used the Gaussian g-binomial coefficients 



-I 9 



(g; g)fc(g; (i)n-k 



involving the g-shifted factorial 

n 

q)n = - -2?""')' ^ e C, n = 0, 1, 
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oo. 



k=l 



We refer to [13] for information about this notation and g-series. We have followed 
the normalization used in Szego [19], where Sq = l/^/q. The Stieltjes-Wigert 
moment problem has been extensively studied in plj using a slightly different 
normalization. 

Lemma 5.1. The double sum 

oo oo 

bj,nbk,nSk+l 

n=0 fc=0 

is absolutely convergent for each j,l > and 



oo 



Moreover, AH = HA = I. 
Proof. We find 



hence for j > k 



I I ^ y 



^— ^ I (T /7 1 in- n\ • 7 , 



< 



-k+p 

{q;q)j-k{q;q)k ^ (g; g)p 
{q'^q)j{q'^q)kf^Q {q;q)p{q;q)oo iq;q)j{q]q)k{q]q)'L'' 



where we have used the g-binomial theorem 



71=0 



with a = = q. By symmetry the estimate holds for all pairs j, k. Since 
Sk+i = it is clear that the double sum is absolutely convergent. 

By Lemma 11.21 we then have 

I = B{B*n) = {I3l3*)n = AH, 

and we clearly have HA = AH. □ 
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From (150!) we get 



P„(0;g) = (-i: 



q2 + A 



hence by (jST 



The matrix /C = (Kj,fc) defined in (fTTl) is given by 



(52) 



(53) 



j+fc 



V(9; g)j(g; q)k ^ 



J 




'k 


.p. 




P. 



hence 



oo 
fc=0 



k=0 ^^'^ p=0 



(54) 



(55) 



in accordance with [1], which also contains other expressions for po- From 
([53]), ([H]) with zo = and [H Theorem 1.2] we get 

1/po < hm Aat < — . 

From the general theory we know that the Stieltjes-Wigert moment sequence has 
an N-extremal solution z/q, which has the mass c = (g; q)ool \^ (=the reciprocal 
of the value in ( |53ll ) at 0. It is a discrete measure concentrated at the zeros of 
the entire function 



D{z) = zY,Pn{^;q)Pn{z- q). 



n=0 



It is also known by a result of Stieltjes, that the measure jl = uq — csq is de- 
terminate, cf. e.g. [El Theorem 7]. The moment sequence (s„) of jl equals the 
Stieltjes-Wigert moment sequence except for the zeroth moment, i.e. 



1 - (?;?)oo)/v^ ifn = 
{n+i)V2 if n > 1, 



q 



and similarly the corresponding Hankel matrices Ti and H. differ only at the entry 
(0,0). 

We shall prove 

Theorem 5.2. The smallest eigenvalue Xn corresponding to the measure jl tends 
exponentially to zero in the sense that there exists a constant A> Q such that 
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The proof of Theorem 15.21 depends on the quite remarkable fact that it is pos- 
sible to find an explicit formula for the corresponding orthonormal polynomials 
which will be denoted Pn{x;q). It is a classical fact, cf. p?, p. 3], that the or- 
thonormal polynomials (P„) corresponding to a moment sequence are given 
by the formula 



/ 



Pn(x) 



n 



det 



So Si • • • 



1 X 



S2n-l 



(56) 



where D„ = det(7i„). In this way Wigert calculated the polynomials Pn{x;q), 
and we shall follow the same procedure for Pn{x; q). Writing 



P„(x; q) = y^&fc.nJ^^, 



A:=0 



we have 

Theorem 5.3. For < k < n 

bk,n = Cn{ — 1 , 



n 

k 

L J g 



where 



C„ 



[-irq 



n I 1 
2 """4 



I.e. 



V(i-(g";g)oo)(l-(g'^+^g)oo)' 

where bk,n denote the coefficients of Pn{x;q). Moreover, 

D„ = D„,(l-(g"+i;g)oo), 
where Dn = det 7i„, D„ = det Tin. 
Proof. We first recall the Vandermonde determinant 



. . . = det 



/ 1 1 

Xi X2 



1 ^n~l n-1 
\Xl X2 



1 \ 



Xtj 



rf.n—1 I 



n I 

l<2<j<n 



(57) 



(5^ 



(59) 



(60) 



(61) 



(62) 
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For an n X n-matrix {aj^k), j,k = 1 . . .n with non-zero elements in the first row 
and column we have 

det(aj-fc) = TTflj,! I TT«i,fc det( ^^'^ ), 

and if aj^k = j = k = 1, . . . ,n + 1, where < g < 1, we get in 

particular 

= det{q-^^^'~'^"/') = (^nV^'/'j det(g-(^-i)('=-i)+i/2), (63) 

hence using Sn = YTj=ij'^ — + 1)(2^ + l)/6 

Dn = g-^"+^+("+^)/V„+i(l, q~\ g-"). (64) 

By (162]) we get 

n n ^ n 

K+i(i,g-\...,g-")=n n -(i-?n=n^~^""^^"^^^'^^'(^'^)"- 



i=0 j=i+l ^ i=0 



and after some reduction 



K+i(l, . . . , g-") = q-'- \{{q. q)r (65) 



We denote by A^+i^p+i respectively A'+i.p+i the cofactor of the entry (r + 1) 
of the Hankel matrix 7i„ = (g^{J+*^^i) Z^) respectively TY^,, where r, p = 0, 1, . . . , n. 
When r = or p = we clearly have A^+i^p+i = A^+i^p+i. For < p < n we get 

An+l,p+l = (-l)"-^det (^g-0-+^-l)V2 I k = l,'.:,n + l\k^r>+l) 
n n+1 



A: = l, ...,n + l;fc7^p+l 



j=l fc=i 



However, 

K+i(l,g-\...,g-") 

p-l n 

= v;(l, q-\ q-^'-'\ g-(^+'\ . . . , g-") n(^"' - ^"'■) n (^"' - ^" 



j=o i=p+i 



V;(l, g-\ . . . , g-(^-^), g-(^+^), . . . , g-")(g; g),(g; g)„_,g-(" +^ 
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so we finally get 

= (-l)"g("+i)(»+V2) ^ ^ , 0<p<n. (66) 

It can be verified that this formula also holds for p = and p = n. 

Using fl56|) it is now easy to verify formula fl50|) for the Stieltjes-Wigert poly- 
nomials P„(a;; q). 

Expanding after the first column we get 

Dn = Dn-cAi^i, c= (g;g)oo/V^> 
and a calculation as above leads to 

which gives (l6Ti) . Moreover, for < p < n we find 

in+i,p+i = An+i,p+i - c(-i)"-pdet (g-o+'^+i)' I fc'==i'':;;n?Mp) , 

and the last determinant can be calculated to be 

-n'^-{n-l)/2+p{p+l/2) 



{q]q)n-p{q]q)p-i^ 
This leads to 

^n+i,P+i ^ ^n+i,p+i 1 - (1 - g)oo 
Dn i-(g"+^;g)oc ' ^ ^ 

It can be verified that this formula also holds for p = because of fl^Tl) . and it is 
now easy to establish (160|1 . □ 

Proof of Theorem 15.21 By Lemma 12.11 we get 

I ..p.. ^^-2 (g;9k(i-(9^+^g)oo)(i-(9^;g)oo) 

Aiv < l^ivlU; gjj ^7v+i/2 • 

From the power series expansion of the entire function {z; q)ao we have 

z 

1 - g 

hence 

iV 

l-(g'^;g)oo~T^, iV-.oo, (68) 
1 - g 



and therefore 



(P^(0;g))-^~^^g^-^^^ iV -> 00, 



which proves the statement of the theorem. 

Remark 5.4. The measure jl is determinate of index 0, cf. |6], so by Corollary 
2.1 in [7J we know that the next smallest eigenvalue Xjq^i of H.^ is bounded below. 
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